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The  theme  of  this  research  has  been  to  develop  effective  strategies  for  the  design  and  analysis  of  flight 
control  systems  for  tailless  fighter  aircraft.  In  our  quest  to  understand  the  fundamental  issues  in  the 
aggressive  maneuvering  of  high  performance  fighter  aircraft,  we  have  developed  tools  and  techniques, 
and  the  appropriate  supporting  theoretical  results,  for  the  study  and  analysis  of  high  performance  ma¬ 
neuvering  systems. 

We  have  developed  a  simplified  flight  dynamics  model,  the  coordinated  flight  vehicle  (CFV).  Exploiting 
the  geometric  structure  of  the  CFV,  one  may  easily  explore  the  space  of  aggressive  flight  maneuvers.  We 
have  developed  an  optimization  based,  trajectory  morphing  technique  by  which  CFV  maneuvers  may 
be  used  to  parametrize  the  achieveable  flight  maneuvers  (i.e.,  maneuvers  that  satisfy  given  high  fidelity 
aircraft  model  dynamics).  We  have  also  developed  theory  and  algorithms  to  allow  these  model  based 
optimization  techniques  to  used  online,  in  a  receding  horizon  fashion.  As  a  theoretical  bonus,  we  have 
obtained  new  results  on  the  structure  of  the  trajectory  manifold  of  a  nonlinear  system. 


Aggressive  Flight  Maneuvers 

It  is  important  to  evaluate  the  performance  of  the  overall  flight  control  system  of  tailless  fighters  over  a 
class  of  aggressive  maneuvers  representative  of  fighter  combat  [8].  Although  a  suite  of  maneuvers  could  be 
supplied  by  recording  the  actions  of  a  pilot  in  simulated  flight,  the  level  of  effort  required  to  equip  a  pilot- 
in-the-loop  simulator  with  appropriate  dynamics  is  difficult  to  justify  in  the  early  stages  of  flight  control 
systems  development  and  evaluation.  Moreover,  it  is  desirable  to  have  methods  for  exploring  the  space  of 
maneuvers  to  identify  potential  shortcomings  in  the  aircraft  or  flight  control  system  design. 

To  this  end,  we  have  developed  a  strategy  for  the  exploration  of  aggressive  flight  maneuvers.  Our  approach 
makes  use  of  the  fact  that  most  aircraft  (and  pilots)  are  designed  to  fly  in  a  coordinated  manner  (i.e.,  without 
sideslip).  Explicitly  enforcing  this  as  a  constraint,  we  have  obtained  a  simplified  flight  dynamics  model,  the 
coordinated  flight  vehicle  [1],  possessing  many  interesting  and  useful  features.  First,  the  coordinated  flight 
vehicle  highlights  the  rich  geometry  of  flight  (e.g.,  the  Hopf  fibration  of  50(3)  pops  up).  Second,  the  CFV 
is  a  nice  control  system  in  that  it  is  differentially  flat  and  path  invertible.  Finally,  the  coordinated  flight 
vehicle  captures  the  global  nonlinear  dynamics  of  maneuvering  flight.  It  is  precisely  these  dynamics  that  a 
pilot  strives  to  manage. 

Using  the  coordinated  flight  vehicle,  we  have  developed  some  prototype  tools  for  the  exploration  and 
visualization  of  flight  maneuvers.  These  tools  utilize  the  OpenGL  graphics  library  (as  a  C/C++  library  under 
Windows  95/NT)  to  allow  3D  dynamic  visualization  (via  animation)  of  an  aircraft  in  maneuvering  flight. 
The  possibilities  are  quite  vast  and  we  are  working  to  develop  an  intuitive  interface  for  control  engineers. 
One  of  the  advantages  of  using  the  coordinated  flight  vehicle  derives  from  its  simplicity — interactive  tools 
can  be  run  on  rather  modest  computer  systems  (e.g.,  PCs  without  graphics  accelerator  hardware). 

One  of  the  tools  we  have  developed  allows  one  to  fly  the  aircraft  by  merely  pointing  in  the  direction  of 
desired  motion.  The  aircraft  rolls  and  climbs  as  needed  to  execute  the  desired  maneuver.  This  is  accomplished 
by  exploiting  the  path  following  properties  of  the  coordinated  flight  vehicle.  With  an  appropriate  interface, 
tools  such  as  this  can  be  used  to  specify  (or,  more  accurately,  develop)  an  aggressive  maneuver  for  evaluation 
on  a  tailless  fighter. 

It  is  interesting  to  speculate  on  the  many  other  potential  uses  of  the  coordinated  flight  vehicle.  Indeed, 
one  can  imagine  using  the  special  geometric  structure  to  guide  the  development  of  autonomous  and  semi- 
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autonomous  flight  management  systems  for  managing  the  global  nonlinear  dynamics  of  flight  vehicles.  In 
that  case,  one  might  truly  develop  a  nonlinear  pilot. 

Now,  the  maneuvers  (or  trajectories)  that  are  generated  using  the  coordinated  flight  vehicle  are  necessarily 
somewhat  idealized.  The  side  slip  angle  f3  is  kept  identically  zero  over  the  entire  trajectory.  Moreover, 
other  important  effects  have  been  intentionally  removed  from  consideration.  For  instance,  the  coordinated 
flight  vehicle  does  not  account  for  the  small  normal  force  that  accompanies  the  generation  of  a  pitching 
moment,  eliminating  the  nonminimum  phase  component  of  the  normal  acceleration  response.  Most  of  the 
low  level  detail  of  actuation  will  also  be  suppressed.  On  the  other  hand,  for  proper  control  system  analysis 
of  maneuvering  aircraft,  one  should  not  shrink  away  from  the  real  world  aspects  of  a  complicated  flight 
dynamics  model,  regardless  of  how  ugly  the  details  may  be. 


Trajectory  Morphing  for  Nonlinear  Systems 

As  noted  above,  we  are  interested  in  finding  trajectories  of  the  complicated  truth  model  even  though  our 
intuition  about  the  nature  of  maneuvering  may  come  from  a  considerably  simplified  system  possessing  nice 
properties.  To  this  end,  we  have  developed  a  homotopy  approach  [3]  for  morphing  a  trajectory  of  a  one  system 
(a  simple  system)  into  a  trajectory  of  another  system  (the  real  system).  These  ideas  will  be  particularly 
useful  in  finding  detailed  trajectories  (including  actuator  dynamics)  for  complicated  systems  such  as  tailless 
fighter  aircraft. 

A  simple  (academic)  example  will  indicate  the  key  ideas  in  morphing.  The  dynamics  of  an  inverted 
pendulum  on  a  cart  are  given  by 

xT  =  f(x,u)T  :=  (x2,sinxi  +  cosx2  u,x4,w)T  . 

where  the  state  vector  is  xT  —  (0, 0,y,y)T  and  the  control  input  u  is  the  cart  acceleration  y.  Although  / 
looks  quite  simple,  it  is  not  easy  to  specify  the  trajectories  of  this  nonlinear  pendulum  system. 

In  contrast,  it  is  very  easy  to  specify,  in  closed  form,  all  of  the  trajectories  of  the  linear  pendulum  system , 

xT  =  f0(x,u)T  :=  (x2,xi  +  u,X4,u)t  , 

given  by  (Jacobian)  linearization  about  the  inverted  equilibrium  position.  Indeed,  since  time  invariant  linear 
systems  are  differentially  flat ,  each  trajectory  of  /0  can  be  specified  by  choosing  a  CA  (or  higher)  function. 
The  flat  output  for  the  linear  pendulum  system  is  y  —  h(z)  =  z$  —  z\.  Note  that  the  nonlinear  pendulum 
system  is  not  differentially  flat. 

Now,  the  general  character  of  trajectories  of  these  two  systems  is  quite  similar  provided  that  the  pendulum 
does  not  depart  too  far  from  the  inverted  position.  Indeed,  we  know  from  linearizing  considerations  that 
for  very  small  pendulum  angles,  the  nonlinear  pendulum  system  will  be  able  to  track  the  trajectories  of  the 
linear  pendulum  system  very  well.  Approximate  tracking  only  works  when  the  desired  trajectory  is  not  too 
aggressive.  There  is  not  much  a  regulator  can  do  when  it  is  fed  the  contradictory  trajectory  information 
derived  from  an  (intentionally)  incorrect  model. 

The  important  point  here  is  that,  even  when  approximate  tracking  is  no  longer  viable,  one  may  be  able  to 
find  trajectories  of  the  more  complicated  system  that  are  much  like  those  of  the  simple  system — provided  that 
essential  couplings  are  similar.  We  exploit  this  idea  and  use  trajectories  of  a  simple  system  to  parametrize 
the  trajectories  of  a  more  complicated  system. 

We  have  taken  the  following  approach.  Given  a  trajectory  ( xo(- ),  ito(*))  of  the  simple  system  /o,  we  seek 
the  trajectory  (xi(*),ui(*))  of  /i  =  /  that  is  nearest  (#o(*)>  uo(*))  h*  a  least  squares  sense.  We  solve  this 
problem  by  tracing  out  a  homotopy  path  in  the  space  of  curves.  In  particular,  (:rx(-),  «*(•))>  0  <  A  <  1,  is 
the  solution  of 

i  J  ||(x(r)  -  x0(r),  u(r)  -  u0(T))||2dr 
Ht)  =  t  6  [0,T] 


minimize 
subject  to 
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where  f\  =  (  1  —  A)/o  +  \  f  defines  a  homotopy  connecting  the  systems.  The  solution  for  A  =  0  is  trivial. 
Under  reasonable  conditions  (involving  differentiation  wrt  A),  the  homotopy  path  in  trajectory  space  will 
be  nice  guaranteeing  a  solution  for  A  =  1.  Moreover,  since  (zaQ? ua(-))  depends  on  A  in  a  smooth  manner, 
the  problem  of  a  good  initial  condition  for  the  optimization  problem  all  but  vanishes.  Figure  1  shows  the 


Figure  1:  Morphed  position  and  angle  trajectories  (A  =  0,0.3, 0.5, 0.7, 1.0). 

morphed  position  and  angle  trajectories  for  the  cart  and  pendulum  system(s)  executing  a  sidestep  maneuver. 
Note  the  angle  trajectories  for  the  linear  and  nonlinear  pendulum  systems  differ  significantly.  This  difference 
is  large  enough  to  make  approximate  tracking  unsuccessful. 


The  Trajectory  Manifold  of  a  Nonlinear  Control  System 

While  developing  trajectory  morphing  techiques,  it  became  clear  that  we  really  didn’t  know  enough  about 
what  the  set  of  trajectories  looked  like.  After  all,  if  we  are  going  to  use  smooth  optimization  techniques, 
wouldn’t  it  be  nice  if  the  domain  of  interest  and  the  functions  involved  possess  some  desirable  smoothness 
properties?  Searching  the  literature,  wre  w'ere  unable  to  find  results  that  set  our  minds  to  rest. 

The  object  of  interest  is  the  set  of  bounded,  exponentially  stabilizable  trajectories  of  x  —  f(x,  u)  which 
we  shall  denote  by  T.  Suppose  that  =  (ao,/^)  €  T  (so  that  doM  =  /(<*<)(*), /* o(*)))  and  that  /  6  Cr, 
r  >  1.  Then,  there  is  a  time- varying  linear  feedback  K  that  stabilizes  fo-  Moreover,  the  projection  operator 
V  :  (a,  p)  »->  (x,  u)  defined  by 

x(t)  =  f(x(t),u(t)),  x(0)  =  a(0) 

u(t)  —  fi(t)  +  K(t)[a(t)  —  x(t)] 

for  t  >  0,  maps  an  L ^  neighborhood  of  £0  onto  the  trajectory  manifold  T.  Researchers  have  used  approximate 
trajectory  tracking  controllers  of  this  sort  for  years.  It  is  a  very  practical  tool. 

The  projection  operator  V  also  turn  out  to  be  a  very  useful  theoretical  tool  [2].  We  have  shown  that 

•  V  is  CT  (wrt  to  Lqo  norm).  Explicit  formulas  for  the  derivatives  are  provided. 

•  T  is  a  CT  Banach  Manifold. 

•  T(T  is  the  set  of  bounded  trajectories  of  the  linearization  of  /  around  £. 

While  results  of  this  sort  are  expected,  there  are  some  surprises.  For  example,  V  is  not  differentiable  in 
L2!  With  a  clearer  understanding  in  hand,  we  have  been  able  to  develop  improved  algorithms  for  local 
optimization  on  the  trajectory  manifold,  including  morphing  techniques. 
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Receding  Horizon  Control  of  Nonlinear  Systems 

The  techniques  developed  above  indicate  that  optimization  can  play  an  important  role  in  the  planning 
and  control  of  high  performance  maneuvering.  For  instance,  once  the  broad  strokes  of  a  maneuver  have 
been  outlined  using  a  coordinated  flight  vehicle  model,  important  details  of  the  low  level  trajectory  for 
tracking  can  be  obtained  using  trajectory  morphing.  To  be  responsive  to  environmental  factors  such  as 
newly  discovered  threats,  it  is  critical  that  the  optimizations  be  conducted  online,  in  a  receding  horizon 
fashion.  Although  online  optimization  has  been  used  extensively  in  the  model  predictive  control  of  chemical 
processes  (relatively  slow  linear  dynamics  with  important  constraints),  its  applicability  to  highly  unstable 
systems  with  fast  nonlinear  dynamics  remained  a  question. 

It  is  well  known  that  unconstrained  infinite  horizon  optimal  control  may  be  used  to  construct  a  stabilizing 
controller  for  a  nonlinear  system.  We  have  shown  [6,  5,  7,  4]  that  similar  stabilization  results  may  be 
achieved  using  unconstrained  finite  horizon  optimal  control.  The  key  idea  is  to  approximate  the  tail  of 
the  infinite  horizon  cost-to-go  using,  as  terminal  cost,  an  appropriate  control  Lyapunov  function  (CLF). 
Roughly  speaking  the  terminal  CLF  should  provide  an  (incremental)  upper  bound  on  the  cost.  In  this 
fashion,  important  stability  characteristics  may  be  retained  without  the  use  of  terminal  constraints  such  as 
those  employed  by  a  number  of  other  researchers.  The  absence  of  constraints  allows  a  significant  speedup 
in  computation.  Furthermore,  it  is  shown  that  in  order  to  guarantee  stability,  it  suffices  to  satisfy  an 
improvement  property,  thereby  relaxing  the  requirement  that  truly  optimal  trajectories  be  found.  We  provide 
a  complete  analysis  of  the  stability  and  region  of  attraction/operation  properties  of  receding  horizon  control 
strategies  that  utilize  finite  horizon  approximations  in  the  proposed  class.  The  guaranteed  region  of  operation 
contains  that  of  the  CLF  controller  and  may  be  made  as  large  as  desired  by  increasing  the  optimization 
horizon  (restricted,  of  course,  to  the  infinite  horizon  domain).  Moreover,  it  is  easily  seen  that  both  CLF  and 
infinite  horizon  optimal  control  approaches  are  limiting  cases  of  our  receding  horizon  strategy. 


Figure  2:  (a)  The  receding  horizon  operating  region  Tj  for  horizon  length  T  =  0.3  and  radius  r  =  rv  =  6.34 
together  with  CLF  operating  region  Hrv.  Also  depicted  are  the  trajectories  x^(-;x)  for  x  on  the  boundary 
of  rJT.  (b)  Receding  horizon  7£(0.3, 0.05)  and  CLF  controller  (dashed)  trajectories. 

The  key  results  of  this  work  have  been  illustrated  on  a  nontrivial  nonlinear  control  system,  the  inverted 
pendulum  on  a  cart,  where  significant  improvements  in  guaranteed  region  of  operation  and  cost  are  noted. 
Figure  2  depicts  the  regions  of  operation  for  the  receding  horizon  strategy  and  a  control  Lyapunov  function 
derived  control  law.  By  comparing  the  closed  loop  trajectories,  we  can  see  that  the  effect  of  the  online 
optimization  is  to  make  better  use  of  system  nonlinearities  providing  significant  performance  improvements. 
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